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ON THE ELASTIC POTENTIAL OF CRYSTALS. 

By William E. Story. 

The theory of elasticity in its applications to homogeneous isotropic 
bodies has been more or less completely developed by Poisson, Cauchy, 
Lame and Clebsch, while Riemann's " Differential-gleichungen " contains an 
admirable abstract of it. 

On a particle of a solid body may act three kinds of forces : 

1) External forces, attractions and repulsions by external masses; 

2) Surface-forces, pressures or tensions applied directly to elements of 
the surface ; 

3) Elastic forces, molecular forces due to the influence of neighboring 
particles of the body itself. 

Elastic forces are produced by unequal displacements of the particles, 
causing a change in their relative positions, and hence a change in the mole- 
cular forces acting between them. Surface-forces affect the internal particles 
of the body only indirectly, by causing a displacement of the particles of the 
surface, giving rise to elastic forces which spread through the interior. It is 
the elastic forces alone that I shall consider in this paper. These may be 
considered as forces tending to prevent any change in the relative positions of 
the particles on one side of a section (plane or curved) of the body with 
respect to those on the other side ; and, because the mutual action of the par- 
ticles becomes inappreciable when the distance between them exceeds a very 
small limit, a particle on one side of the section will be affected only by those 
particles on the other side which lie very near the section and within a small 
distance of the particle in question. So that the whole mutual effect of the 
particles on the two sides of the section may be treated as the sum of the 
mutual effects of particles on the two sides of the very small portions (which 
may be considered plane-elements) of the section ; and it is only necessary to 
consider the elastic forces on such plane-elements, which vary from element to 
element of the section, but, in the immediate vicinity of any given point of the 
section, are proportional to the plane-elements on which they act. It is there- 
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178 Story, On the Elastic Potential of Crystals. 

fore convenient to call the components of the elastic force acting on any plane- 
element, the ratios of the actual components to the area of the element. The 
components of the elastic force acting on any plane-element passing through 
a given point may be expressed in terms of the components of the elastic 
forces acting on any three mutually perpendicular plane-elements through the 
point, e. g. three elements parallel to the planes of a rectangular system of 
coordinates. 

Let x, y, z be the coordinates of any point of the body referred to a fixed 
rectangular system, and let three plane-elements of areas e x , e y , s z , whose nor- 
mals are parallel to the axes of x, y, z respectively, pass through this point ; 
then in notation of Cauchy the components of the elastic forces acting on 
e x , s p , e z in the directions of the three axes will be X x , Y x , .Z x ; X y , Y y , Z y ; 
X z , Y z , Z z : in each of which the large letter denotes the direction of the com- 
ponent, and the suffix the direction of the normal to the plane-element on 
which it acts. Let, further, the density at the point x, y, z be p. 

If x, y, z is a point in the interior of the body, let the components of 
the resultant of the external forces acting on a very small mass m at this 
point be ml, mY, mZ; then the conditions for the equilibrium of this small 
mass are 

v bX x . hX y bX s A 
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If, however, x, y, z is, a point of the surface, let a, (3, y be the direction- 
cosines of the inner normal to the surface, and eB, sH, sZ the components of 
the surface-force on an element e of the surface ; then the conditions for the 
equilibrium of an infinitesimal solid element immediately under the surface- 
element e are 

(2) ] H+aY, + pY, + yY. = 0, 

( Z+aZ. + pZ 9 + yZ, = 0, 

where the surface-force is considered positive if a pressure and negative if a 
tension. 
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It is well known that, if u, v, w are the displacements, in the directions of 
the axes, of the point x, y, z, and if 
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then there exists a function <5 of x x , y y , z z , y z , z x , x y , such that 

(4) 
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This function 4>, which may be called the elastic potential, is homogeneous 
of the second degree in x x , y y , z 2 , y z , z x , x y , and contains therefore in general 21 
terms, whose coefficients are characteristic for the body under consideration, 
depending only on its structure in different directions. For a homogeneous 
body, i. e. one, every particle of which is similarly surrounded by particles, 
the coefficients of <£ are constant, but for a non-homogeneous body they differ 
from point to point, i. e, are functions of the coordinates x, y, z. For a homo- 
geneous isotropic body the form of <J> has been determined to be 

(5) $ = i A, {x x + y y + z z ) 2 + {i (x 2 x + y\ + z\) + ~ (i (yl + z\ + a$ > 

where /I and ^ are constants. I now propose to determine its form for homo- 
neous crystalline bodies. 

Each crystalline system has a form of symmetry with respect to three 
mutually perpendicular axes which is not altered by a transformation to any 
one of certain other systems of mutually perpendicular axes having the same 
origin. We may then take one of these systems of axes as the basis of a 
system of coordinates x, y, z, and any other of them as the basis of a system 
of coordinates x', y 1 , z'. The external form of symmetry of a crystal must, it 
would seem, be also a form of symmetry of arrangement of its particles ; if 
this be so, a transformation of coordinates from x, y, z to x 1 , y", z' must leave 
the coefficients unchanged. This will be the basis of the following determina- 
tions of the relations between the coefficients, on which relations the form of 
4> depends. 

The transformation from x, y, z to x', y, z' can be effected by a combina- 
tion of two simple kinds of transformations, viz. by inversions of axes and 
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rotations about axes through certain angles. There are six such transforma- 
tions, which I will designate by the first six letters of the alphabet, as follows : 

a) inversion of the axis of x , 

x 1 = — x , y" = y , z' = z , u = — vl , v = tf , w = w' , 

x x = *£/ , y y = y y f , z 2 = z 2 / , y a = y z t , £„ = z x * , # y = ^^ ; 

b) inversion of the axis of y , 

x 1 = x , y" = — y , z' = z , u = u' , v = — v' , w = w' , 

# = .£3;' , y y = Ifyl , Zj, = Zgf , J/ 2 = y z l , £,. = 2^/ , Xy = Xyl ) 

c) inversion of the axis of z , 

x' = x , y 1 = y , z' = — ■ z , u = u' , v = v' , w = — w' , 

x x =■ x % i , y y = y y / , 2-j, = 2 y , y,j = y .j> , £,. = z % i , ^ •= * ,j f 

d) rotation about the axis of x through the angle a , 

x 1 = x , y" = y cos a -\- z sin a , z' = — y sin a + z cos a , 

u = u' , v = v' cos a — w' sin a , w = i/ sin a + w ' cos a > 

x x = ^ , y y =y' y , cos 2 a — yy sin a cos a + *v sin 2 a , 

3 2 . = y"y, sin 2 a + yV si* 1 a cos a + z' z , cos 2 a , 

y. = (l/v — z '*) sin 2a + yV cos 2a , 

^ = afy, sin a + z' xJ cos a, x P = af^ cos a — z'^ sin a ; 

e) rotation about the axis of y through the angle /3, by formulae found from 
those of d) by changing x to y, y to z, z to #, and a to /3 ; 

/') rotation about the axis of z through the angle y, by formulae found from d) 
by changing x to z, y to x, z to y, and a to y. 
The general form of <J> is 

$ = a n + # 2 + a i2 y y + « 33 z 2 + 2ff 23 y 2 ,z, + 2a 13 x x z z + 2a 12 ^yj, 
+ 2a u x x y, + 2a 15 a?A + 2a 16 ^^ + 2a 2i y y y z + 2a 9 jy y z x 
+ 2a w yyXy + 2a 34 z,y, + 2^2^^ + 2a 36 z^ + a u yl 
+ a ss z 2 x + a 66 o; 2 + 2« 56 z^ + 2a m x y y z + 2a45y 2 2;, . 

That the transformation a) shall leave $ unaltered in form it is necessary 
and sufficient that 

( / ) @>\5 — #25 == ^35 == ^45 := *56 := ^16 == #26 — #36 == #46 == "• 

That the transformation b) shall leave the form of <f> unaltered it is necessary 
and sufficient that 

(o) <?14 — <?24 -— #34 := #45 := #16 = — #26 — #36 — ■ #56 -— "• 

That the transformation e) shall leave the form of <j> unaltered it is necessary 
and sufficient that 

(9) a u ~ a u — a u = a 46 = a w = a w =z a w = « 56 = 0. 



(6) 



(10) 
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That the transformation d) shall leave the form of <E> unaltered it is necessary 

and sufficient that 

(a 12 — a 13 ) sin a = , (# 55 — # 66 ) sin a = , (# 22 — # 33 ) sin a = , 
a 14t sin a = , # 5G sin a = , a ls (1 — cos a) = , a m (1 — cos a) = , 

(# 22 + #33 2«23 4«44) SHI a COS « = , 

(# 24 + #34) sin a = , (# i4 — # 34 ) sin a cos a = , 

(# 25 + # 35 ) (1 — cos a) = , (#, 6 + # 36 ) (1 — cos a) = , 

[(#25 — %>) + (#26 — #35)] sin a cos a + ( cos a — cos 2a — sin a) $45 

+ (cos a — cos 2a + sin a) # 46 = , 

[(#25 — #35) — (#26 — <%)] sin a cos a + (cos a — cos 2a + sin a) a m 

— (cos a — cos 2a — sin a) # 46 = , 
(a 25 — %,) (cos a — cos 2a — sin a) + (« 26 — #35) (cos a — cos 2a + sin a) 

— 4 (a & + #4e) sin a cos a = , 
(a 25 — <%) (cos a — cos 2a + sin a) — (# 26 — # 36 ) (cos a — cos 2a — sin a) 

— 4 (# 45 — #«) sin a cos a = . 

The conditions (10) are, of course, satisfied by a = 0. They will also be 
satisfied by a multiple of 90° when certain relations hold between the 
coefficients, as follows : 
if a = 90° or — 90°, 

# 14 ~ a 1B ~ ff 16 ~ $2 5 rr # 26 = #35 — #36 — #45 — #46 — #56 — U > 



C 11 ) , 

#12 #13 ) #22 #33 ) #55 #66 > #24 #34 5 

if a = 180°, 

(12) #ig — #25 — #35 = #45 — #16 == #26 — #36 — #46 = " ) 

which are identical with (7), so that, if an axis can be inverted without pro- 
ducing any effect on <3>, a rotation about that axis through an angle of 180° 
will also produce no effect, and vice versa, as is also evident from the fact that 
all three axes may be inverted without producing any effect on the most gene- 
ral form (6) of <J>. 

If a has neither of the above values, 

#14 — #24 — #34 — #15 — #16 — #56 — " ) 



(13) 

1 # 12 1= #13 , #22 — #33 — #23 ~f~ 2#44 , #55 — #66 ) #25 = #35 5 #26 

and either 

(13i) #25 = #35 — #45 — #26 — #36 — #46 — O , 

or 

(13 2 ) 4 cos 6 a + 4 cos 5 a — 2 cos 4 a — 2 cos 3 a — 4 cos 2 a — 4 cos a — 1=0 
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(whose only real roots between — 1 and + 1 are cos a = — 0.607 + and 
cos a ■=. — 0.389 — ), i e. if a rotation about an axis through any other angle 
than 90°, 180°, 270° (or — 90°), cos [ - 1] (— 0.607 +) or cos^ 1 ^— 0.389 — ) 
does not affect the form of <&, a rotation about this axis through any angle 
whatever will not affect it, i. e. there exists perfect symmetry about this angle. 

The transformations e) and /) give conditions similar to (10), (11), (12) 
and (13), found from them by a cyclic interchange in each set of suffixes 1, 
2, 3 and 4, 5, 6 and in a, (3, y. 

Each of the crystalline systems admits of a characteristic combination of 
the above transformations. 

The triclinic system admits only of an inversion of all the axes at the 
same time, which leaves the most general form of $ unchanged without any 
conditions between the coefficients, i. e. (6) is the general form of <I> for tri- 
clinic crystals. 

The monoclinic system, in which the axis of x is taken in that axis of the 
crystal which is perpendicular to the plane of the two others, allows only the 
transformation a) and d) in which a — 180°, i. e. its conditions are (7) and 
(12), which are identical, and 

( <3> = a n xl + a 2 . 2 f y + a 3S zl + 2a 23 y y z z + 2a n x x z z 

(14) J + 2a 12 x x y y + 2a u x x y z + 2a u y y y z + 2a u z z y z + a u yl 

( + a BS zl + a^Xy + 2a m z x x v . 

The trimetric system, in which the axes of x, y and z coincide with the 
axes of the crystal, admits of each of the transformations a), b) and c), hence 
for this system the conditions are (7), (8) and (9), and the form of <I> is 
$ = a u xl + a 22 y 2 v + a m zl + 2a^y y z z + 2a n x x z z + 2a n x x y y 



( 15 ) ] - , ./> , 

The dimetric system, whose principal axis is the axis of x and whose 
secondary axes are those of y and z, admits of the transformations a), b), c) 
and d) where a rr 90°, so that the conditions are (7), (8), (9) and (11), and the 
form of $ is 

$ = a n xl + « 22 (y y + zl) + 2a & y y z z + 2a n x x {y y + a.) 



{ + UuVl + «55 {Zl + Xy). 

The monometric system, whose axes are the axes of x, y and z, admits of 
the transformations a), b), c), d), e) and/') for a =. 90°, (3 — 90°, y = 90° ; so 
that the conditions are (7), (8), (9), (11) and, similar to (11), 
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#14 = #16 = #24 — #25 — #26 — '■ #34 — #36 — #56 — #46 — #45 — " ) 

#23 = #12 ) #11 — #33 ) #44 — #66 ) #15 — #35 ) 

#14 — #15 — #24 = #25 — #34 — #35 — #36 — #56 — #46 — #45 — " ) 
#23 — #13 ) #11 — #22 ) #44 == #55 ! #16 — #26 ) 

and the form of <$> is 

(17) $ = #n (4 + $ + «J) + 2a u (&#, + ^ + ^) + ^ ($ + z\ + 4). 

The hexagonal, as well as awy regular prismatic system, other than one 
whose secondary axes include an angle of 90° (because the conditions (13i) 
omitted in case a = cos [_1] ( — 0.607 +) or a = cos [_11 ( — 0.389 — ) are more 
than replaced by (7) ), whose principal axis is that of x, admits of the trans- 
formations a), b), c), d) ; the conditions are then (7), (8), (9), (13) and (13j) 
and <I> has the form 

f * = #n4 + #22 {f y + zl) + 2a 2S y y z z + 2a l2 x x (y y + z t ) 

( +4 K* — a ^ $> + ^55 (*2 + X D ■ 

It is to be remarked that (16) differs from (18) only in the coefficient 
of yl, which is independent of the other coefficients in the former case, but 
dependent on the coefficients of y y + z\ and y y z z in the latter case, so that, as 
far as their elastic relations are concerned, square prismatic crystals seem to 
be the most complex of regular prismatic forms, in as much as they contain 
one more constant than any other. 

From the above determinations of <£ it will be seen that the elastic poten- 
tials of isotropic bodies, of monometric, hexagonal, dimetric, trimetric, mono- 
clinic and triclinic crystals contain respectively 2, 3, 5, 6, 9, 13 and 21 con- 
stants. The forms of <S>, (5), (17), (18), (16), (15,) (14) and (6) substituted 
in (4) give the forms of the components of the elastic force in these differ- 
ent cases of homogeneity, - and these components substituted in (1) and 
(2) give the conditions for equilibrium under a given set of external forces 
X, Y, Z and surface-forces B, S, Z. The conditions for motion, e. g. elastic 
vibrations, will be found from those for equilibrium by the substitution 

(fi/ii ffiv ffiw 

of X — , Y ^— —,, Z — r-r- for X, Y, Z respectively, t denoting time. 

air dr dir 



